Quantum Noise in the Collective Abstraction Reaction A+B 2 ^AB+B 
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We demonstrate theoretically that the collective abstraction reaction A+B2 — > AB+B can be 
realized efficiently with degenerate bosonic or fermionic matter waves. We show that this is domi- 
nated by quantum fluctuations, which are critical in triggering its initial stages with the appearance 
of macroscopic non-classical correlations of the atomic and molecular fields as a result. This study 
opens up a promising new regime of quantum degenerate matter-wave chemistry. 
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The making and probing of ultracold molecular gases 
has attracted much attention in recent years [l[ . Starting 
from an atomic Bose-Einstein condensate, magnetic Fes- 
hbach resonances and/or optical photoassociation 0] can 
be exploited to create not only diatomic molecules but 
also more complex molecules, such as evidenced by the 
recent observations of transient short-lifetime trimcrs CS3 
[H or tetramers CS4 (actually the resonances in inelas- 
tic processes) In another development of relevance 
for the present study, atom-dimer dark states were pro- 
duced through coherent photoassociation [||, a process 
sometimes called superchemistry or quantum degenerate 
chemistry Q. 

So far, matter- wave superchemistry has concentrated 
largely on the coherent combination or decomposition 
reactions |7f between atoms and homonuclear or het- 
eronuclear [8( molecules. In this Letter we extend these 
ideas to another type of elementary chemical reaction, 
the coherent abstraction reaction (or bimolecular reac- 
tive scattering) A+B2~>AB+B. This reaction is an im- 
portant benchmark system extensively studied for many 
years in chemical physics, a particularly noteworthy con- 
tribution being the study by Shapiro and Brumer of the 
coherent control of single-molecular photoassociation or 
bimolecular collisions via interference of reactive path- 
ways [Ij] . One main result of this Letter is to demon- 
strate that the coherent abstraction reaction can be real- 
ized and controlled efficiently in degenerate matter waves 
by exploiting an atom-molecule dark state. An impor- 
tant characteristic of this process is that it is triggered 
by quantum noise, leading to large shot-to-shot fluctu- 
ations that dominate the initial stages of the reaction. 
From a theoretical point of view, this implies that the 
mean-field Gross-Pitaevskii equation is not appropriate 
to describe the early stages of the coherent bimolecular 
reaction, and can be used only once the product reac- 
tant populations become macroscopic. This is in contrast 
with single-molecular combination reactions such as the 
atom-dimer 0, E| or atom-trimer conversion [l~oj ] , and is 
reminiscent of quantum or atom optics situations such as 



the laser, superradiance and matter-wave superradiancc, 
see for instance Refs. [Il, 12, 13 . 

The basic idea in realizing the collective reaction A+ 
B2 -^AB+B in quantum-degenerate gases is to first cre- 
ate highly excited trimers AB2 via an entrance-channel 
atom-dimer Feshbach resonance, and to then photo- 
dissociate them into a closed-channel bound dimer and 
an atom. The use of a dynamical two-photon resonance 
scheme involving an intermediate trimer state permits 
one to exploit the existence of a coherent population trap- 
ping state (CPT) that prevents the trimer population 
from becoming significant throughout the conversion pro- 
cess. Such a generalized atom-molecule dark state does 
not exist in other schemes that involve e.g. an intermedi- 
ate two-species atomic state. Note also that this scheme 
is different from a purely collision-induced reaction 
and from the nondegenerate single-pair dynamics of re- 
active scattering To our knowledge this is the first 
proposal for the quantum control of matter- wave abstrac- 
tion reactions, and as such it represents a promisingnew 
step in developing the field of degenerate chemistry [J, 0] . 
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FIG. 1: (Color online) Schematic illustration of the coherent 
abstraction reaction A+B2 — >AB+B with degenerate matter 
waves. Here the atoms A and B can be bosonic or fermionic. 

Our model system is illustrated by Fig. 1. Denoting 
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the strength of the A+B2^AB2 coupling with detuning 
S by X[, the Rabi frequency of the dissociating laser by 
ft^ and its detuning by A, the dynamics of the system is 
described at the simplest level by the model Hamiltonian 



n 



+ X' 
- ft' 



^](r)^ a (r)^ b2 (r) + h.c. + (A + 6)^ ab (r)^ ab (r) 



(i) 



We consider first the case of a bosonic system, in which 
case the annihilation operators ipi, where the indices 
i, j = a, b, 6 2 , ab, t stand for atoms (A and B), dimers (B 2 
and AB) and trimers, satisfy standard bosonic commu- 
tation relations. The terms proportional to \'i j describe 
s-wave collisions between these species. We remark that 
trimer formation via an atom-dimer resonance is actively 
studied in ongoing experiments, and that the Feshbach- 
resonance-aided photoassociation considered here could 
also be reformulated as a laser frequency modulation 
scheme (j| [l5j . Finally, the Bose-enhanced selectivity of 
dissociation channels in coherent photodissociation of the 
heteronuclear trimers ABC was also studied by Moore 
and Vardi Q. 

In the framework of a standard mean-field approach 
where ipi — > ^/nipi we find readily 
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with Xi,j = n Xi.ji A = X'^fn and ft = ft'y 7 ", and the 
phenomenological decay rate 7 accounts for the loss of 
intermediate trimers. The maintenance of dark or un- 
populated trimer state is important for an efficient con- 
version even with a very short molecular lifetime [HI]. 
For -06(0) = ^ab(O) = 0, Eqs. (2) imply that there is 
no growth in the populations of the atoms B and of the 
dimers AB. This indicates that the mean-field Gross- 
Pitaevskii equations, which are adequate to describe the 
short-time behavior of the familiar coherent atom-dimcr 
or atom-trimer conversion 0, [To| , break down completely 
in studying the onset of this bimolecular reaction. 

A similar situation has been previously encountered in 
a broad range of systems in quantum optics, but also in 



coupled degenerate atomic and molecular systems such as 
the matter-wave superradiance in Bose-condensed atoms 
12, llj]. Following a strategy developed in the study 
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of these systems, we decompose the problem into an ini- 
tial stage dominated by quantum noise followed by a clas- 
sical stage that arises once the product components have 
acquired a macroscopic population. The initial quantum 
stage is treated in a linearized approach whose main pur- 
pose is to establish the statistical pro per ties of the initial 
fields required for the classical stage [13l |. 

To simplify the description of the initial stages we 
note that in the collisionless limit an effective second- 
quantized Hamiltonian can be obtained by adiabatically 
eliminating the intermediate excited state, 



7-Lctf = -{Gc\ b c\c a Cb 2 + h. 



c.) + c , 



(3) 



where c = uJic a c a c\ 2 c b2 



^2C ab c ab c\c b , 



G = (X'il'/S) J dr^rWb{r)Ur)hAr), 

wi = (A' 2 M) j dr<t>* a (r)ct>a(r)<t>t 2 (r)<f>b 2 (r), 

"2 = (ft' 2 /<5) / dr^ ab (r)ct> ab (r)^Ur)Mr), 

and tpi(r,t) = 4>i( r )ci(t) ■ Equation (3) is reminiscent of 
the Hamiltonian describing spin-exchange scattering in a 
two-species two-pseudospin-state Bose condensate [16| . 

For short enough interaction times, the populations of 
the products remain small compared to the total parti- 
cle numbers. In this regime, we can treat the fields tj) a 
and i\)b 2 classically, c a ,b 2 —> \/N a ,b 2 , and then neglect the 
term in Eq. (3) describing only the interactions of the 
modes c a ,b 2 - This amounts to linearizing the dynamics 
of the fields c ab and c b , with the noise source fj(t), 



Ub(t) = fl ab (t) = igci ab (t), 



(4) 



which is familiar from quantum treatments of the optical 
parametric oscillator [lll| and of the molecular dissocia- 
tion (pair production) [171 ] . The noise operators satisfy 

(fHt)f j (t')) = o, (f i (t)f](t')) = g% j 6(t-t'). 



where Q = Gy/ N a Nb 2 and i,j — ab or b here and in 
the following. It is these noise operators that trigger the 
evolution of the system from initial vacuum fluctuations. 

The quantum noise- induced populations of the c g2j b 
mode and their fluctuations correlation are 

Nj = (c]dj) = sinh 2 (gt) ?z N a N b2 G 2 t 2 ; 

(AN ab AN b ) . 2 

Cab =C b = = 1 + smh (gt) > 1, (5) 

V^ab^b 

where ANj = Nj — (Nj}. Equation (5) can also be de- 
rived by solving Eq. (3) to second order in time. The 
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second factorial moment of the modes £j is typical of 

(2) 

chaotic fields, = 2, but they are entangled, with 

2 
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0.035 



(2) 
9ab.b 



,(2)„(2) = sn^r 2 {g t ) + Um\r x {gt) > 0, (6) 



indicative of a violation of the classical Cauchy-Schwartz 
inequality. 

For comparison we comment briefly on the case where 
atoms A are bosonic and atoms B fermionic. We obtain 
similar equations of motion in that case, except that the 
noise operators (—ft, P ab ) are now different. These equa- 
tions can be solved via a Bogoliubov transformation, and 
wc find that as a result of the Fermi statistics the vacuum- 
noise-triggered population is now Nj = sin 2 (Qt) < 1 
(for the zero- momentum mode fl7l]). with the dimer-atom 
pairs correlation becoming C a b = C b = 1 — sin 2 (C?<) < 1. 
The Mandel Q parameters are 



(a) cosh 2 (gt) > 1, 

(b) cos 2 (Gt) < 1, 



(7) 



where (a) is for creating bosonic and (b) fermionic 
matter-wave fields, which exhibit therefore super-Poisson 
or sub-Poisson statistics , respectively. 

We conclude the discussion of the short-time dynamics 
by mentioning that the long-time quantum statistics of 
the AB and B populations, which are significantly influ- 
enced by the initial vacuum fluctuations, can be calcu- 
lated by a positive-P representation technique [l8| and 
other methods [lj|. Rather than adopting such a full 
quantum treatment, we proceed in the following by com- 
bining the mean-field description of Eqs. (2) with stochas- 
tic classical seeds with statistics consistent with the re- 
sults of the short-time linearized quantum theory [l]| . 

We now turn to the long-time reaction dynamics. We 
proceed by numerically computing a large number of tra- 
jectories (typically about 300) from initial classical seeds 
that satisfy the short-time statistics of Eq. (5). For each 
trajectory n, we use Eqs. (2) to calculate the particles 
populations Ni >n (t) where i = A, B2, AB, and B. 

Figure 2 shows the standard deviation ANi(t) of the 
particle populations i=A, B, AB and B2 for 5 = ±3, 
and the insert shows a range ±ANi about their mean, 
Ni(t)±ANi(t) for S — 3 and in the case of bosonic atoms. 
As expected, the small product seeds triggered by the 
initial quantum fluctuations are significantly amplified 
before reaching a stationary value for S = 3. 

An important feature of the coherent abstraction reac- 
tion is that it can be controlled and optimized by exploit- 
ing the existence of a CPT dark state 0, [lj| • Under a 
dynamical two-photon resonance condition, Eqs. (2) ad- 
mit a steady-state CPT solution such that a trimer state 
remains unpopulated at all times 
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FIG. 2: (Color online) Standard deviations of the atom and 
dimer populations (ANj(t) = (AiV J 2 } 1/2 ) for S = 3 or 5 = -3. 
Time is in units of A -1 , and 7 = 1. The other parameters 
are given in the text. The trimer number remains zero at all 
times due to the CPT condition. Inset: fluctuating range of 
the populations Ni(t) ± ANi(t) for 5 = 3. 



where we have applied the steady-state ansatz [lj| V*| = 
l^le^e^ 4 , (6,rib 2 ,ab = 2(9,n) h ,{6,n) a = (0,fi) b and 
7Z = N a (Q)/2Nb 2 (0) is introduced to define the initial 
ratio of the particles numbers. Using dN^ b /dlZ = 0, we 
can find the maximum value 2V* b | max = 1/3 for 1Z = 1/2. 

Figure 3(a) shows the mean particle populations ob- 
tained by averaging over tit = 300 trajectories. In this 
specific example atom A is 87 Rb and atom B is 41 K, A = 
4.718 X lOS" 1 and fi(f) = fi scch(t/r) with Q a /X = 20, 
Ar=20 20]. The collision parameters, in units ofX/n, are 
Xaa = 0.5303, Xbb - 0.3214, X ab = 0.8731, and the oth- 
ers arc taken as 0.0938 [20l | . We note that the scattering 
lengths of the various collisions, especially those involving 
molecular trimers, arc not known at this time. Wc have 
therefore carried out numerical simulations by for a large 
set of plausible collision parameters for the Rb-K, Rb-Na 
or other atomic condensate [3] . We found that the stable 
bimolecular conversion is always possible for appropriate 
values of the external field detuning S. The departure of 
the product populations from the ideal CPT value is due 
to the fact that only an approximate adiabatic condition 
exists for the CPT state [lj| \j<%: j nl (t) w J^L^ < 1, 
with rj = A/0, which becomes increasingly difficult to 
satisfy in the last stages (rj ^> 1). 

Our proposal relies crucially on the capability to avoid 
rapid collisional quenching or the formation of an un- 
stable atom-dimer sample. When energetically allowed, 
collision-induced reactions always occur at some rate, 
and we need to guarantee that the time scale over 
which quantum fluctuations dominate the dynamics is 
short enough, so that the dynamics of the system is not 
collision-dominated. From the condition \Qt\ < 1 we es- 
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FIG. 3: (Color online) Populations of dimers and atoms for 
S = 3 or 5 — —3 (in units of X/n) by averaging 300 trajecto- 
ries. Time is in units of A -1 , and 7 = 1. The trimer number 
remains zero at all times. The line labelled "CPT" shows the 
ideal population of products (dimers AB and atoms B). 



timate the maximum permissible collision-induced rate 
to be of the order of 10 5 s" 1 for |<5| = 3 and fi = 20A. 
A typical low-temperature inelastic collisions rate coeffi- 
" 17 m 3 /s [H, 21|, corresponding to a rate of 
for a sample density of 10 14 /cm 3 . In that 



cient is 10 
about 10 3 s _1 

case the fluctuations-induced dynamics will dominate the 
short-time behavior of the system. 

Finally we remark that by using 87 Rb for the bosonic 
atoms A and 40 K for the fermionic atoms B we can realize 
the conversion of bosonic pairs to fermionic pairs, which 
is likewise described by Eqs. (2) with the substitution 



^2/3 . 



|4/3 



(9) 



where A 3 = (Qir 2 ) 2 ' d /4Mj, j = ab,b, and Mj denotes 
the particle mass, provided that we ignore the s-wave 
collisions of fermionic particles and only consider their 
kinetic energy [22J. In contrast to the purely bosonic 
case, molecule formation is now found to be stable for 
both positive and negative detunings, see Fig. 3(b). 

In conclusion, we have shown that in the collective ab- 
straction reaction A+B2 — ► AB + B, the initial quantum 
fluctuations lead to the strongly correlated creation of 
dimer-atom pairs. This novel feature indicates that in 
the quantum degenerate regime, elementary bimolecu- 
lar reactive scattering are significantly different from the 
familiar single-molecular reaction, and may lead to fasci- 
nating new opportunities in degenerate chemistry, such 
as the possible collective reaction 2A2 — > A3 + A. 



Future work will also study the unique "supcrchem- 
istry" effects of ultra-selectivity or conhnement-induced 



stability in our system 
sis of collisional effects 




In addition, a complete analy- 
141 [23| will need to be considered. 
While experiments along the lines of this analysis promise 
to be challenging, recent progress in quantum degenerate 
chemistry [5] . I a ha and in the control of atom-molecule 
systems |3j, U, [2J, [25J indicates that achieving this goal 
should become possible in the not too distant future. 
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